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Abstract 



| In order to establish links between p-wave pion production in nucleon-nucleon collisions and low energy three-nucleon scattering, 
j an extensive programme of experiments on pion production is currently underway at COSY-ANKE. The final proton pair is 
. . detected at very low excitation energy, leading to an S-wave diproton, denoted here as {pp} s - We now report on measurements of 
^ the differential cross section and analysing power of the pp {pp} s 7r° reaction at 353 MeV. Both observables can be described 
in terms of s- and d-wave pion production and, by using the phase information from elastic pp scattering, unique solutions can 
K ^N be obtained for the corresponding amplitudes. This information is vital for the partial wave decomposition of the corresponding 
-~ pn — > {pp} s 7r~ reaction and hence for the extraction of the p-wave terms. 
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Within the context of chiral perturbation theory, a sig- 
nificant step forward in our understanding of pion physics 
at low energies would be to establish that the same short- 
ranged NN — > NNtt vertex contributes to p-wave pion 
production, to low energy three-nucleon scattering [1—3], 
7<i — > ttNN [4,5] and ird — > jNN [6], as well as in weak 
reactions like tritium beta decay [7-10]. The relevant tran- 
sition amplitude, which connects NN S'-waves in the ini- 
tial and final state with a p-wave pion, contributes to both 
pp — > 7r + d(n + pn) and pn — > ppix~ . However, the extensive 
data for 7r + production is of limited use in this context, be- 
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cause the p-wave amplitudes are completely dominated by 
the 1 D2 initial state, which hinders a reliable extraction of 
the 1 S'o initial state [3] . 

There is a programme at the COSY-ANKE facility of 
the Forschungszentrum Jiilich to perform a complete set 
of measurements on NN — > {pp} s ir at low energy [11]. 
Here the {pp} s denotes a proton-proton system with very 
low excitation energy, E pp . At ANKE we select events with 
E pp < 3 MeV and, under these conditions, the diproton 
is overwhelmingly in the 1 Sq state with antiparallel pro- 
ton spins. This simplifies enormously the spin structure: a 
partial wave analysis for pp — > ppir without the E pp cut 
would require twelve additional P-wave final pp spin-triplet 



states [12-14]. The cut also allows one to extract the full 
information on the production amplitudes without having 
to make measurements of the final proton polarisations. 

Whereas the {pp} s final state is isotriplet, the isosinglet 
np initial state also contributes to pn — > ppir~ . In order 
to isolate this, which contains the amplitudes of interest, 
the isotriplet channel needs to be well understood. As the 
first part of the outlined larger programme, we therefore 
report here on measurements of the cross section and pro- 
ton analysing power in the pp — > {pp} s t° reaction at T p = 
353 MeV. 

For a spin-singlet diproton, the spin structure of 
the pp —¥ {pp} s 7T° or np —> {pp} s ir~ reaction is that of 
i ^ — > + _ . Parity and angular momentum conserva- 
tion require that the initial nucleon-nucleon pair to have 
spin S = 1. The pion orbital angular momentum I and 
the initial nucleon-nucleon isospin I are then linked by 
I + I = odd so that, for the pp — > {pp} s n° reaction, only 
even pion partial waves are allowed. As a consequence, the 
unpolarised cross section for 7r° production, and this times 
the proton analysing power Ay, must be of the form 

da\ k 

dn) = ip + a2 cos2 &7r + ° 4 cos4 ^ H — ^ 

A y (Jj^J = ^ sin 6 n cos e n (6 2 + 64 cos 2 9 n H ) , (2) 

where 9^ is the pion cm. production angle with respect to 
the direction of the polarised proton beam. Here p is the 
incident cm. momentum and k that of the produced pion 
which, at 353 MeV, have values p = 407 MeV/c and k « 
94 MeV/c, where the latter represents an average over the 
3 MeV Epp range. 

The only detailed measurements of the pp — > {pp} s tt° 
differential cross section over the whole angular range were 
carried out with the PROMICE-WASA apparatus at CEL- 
SIUS at a series of energies from 310 to 450 MeV, using 
the same standard 3 MeV cut on E pp [15]. Throughout this 
energy range, significant anisotropies were found in the an- 
gular distributions which were attributed to interferences 
between pion s and d waves. On the other hand, there were 
no corresponding measurements of the proton analysing 
power, which might also be driven by a strong s-d interfer- 
ence. 

We have previously reported measurements of the 
pp — > {pp} s 7T° differential cross section at several energies 
and small angles [16,17]. Since these were carried out using 
the ANKE spectrometer [18] under conditions that were 
similar to the current ones, the description here can be 
quite brief. ANKE is placed at an internal beam station 
of the COSY cooler synchrotron. Fast charged particles, 
resulting from the interaction of the stored transversally 
polarised proton beam with the hydrogen cluster-jet tar- 
get [19] and passing through the analysing magnetic field, 
were recorded in the forward detector (FD) system. The 
FD, which was the only detector used in this experiment, 
includes multiwire proportional chambers for tracking and 



a scintillation counter hodoscope for energy loss and timing 
measurements. 

To start the identification of the pp —t {pp} s it reaction, 
proton pairs were first selected from all the registered two- 
track events using the measured momenta of the both par- 
ticles and the difference in their time-of- flight [20] . The res- 
olution a(Epp) in the diproton excitation energy was better 
than 0.6 MeV, which allowed the E pp < 3 MeV cut to be 
applied reliably. 
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Fig. 1. Two-dimensional distribution of the missing-mass-squarcd 
Mv of the pp — > {pp}sX reaction at 353 MeV versus the diproton 
cm. polar angle (9™ for events with E pp < 3 MeV. 

After selecting the 1 Sq final state, the kinematics of the 
PP —> {pp}sX process could be reconstructed on an event- 
by-event basis to obtain a missing-mass Mx spectrum. A 
two-dimensional distribution of M x versus the cm. polar 
angle of the diproton #™ is presented in Fig. 1. This demon- 
strates the large angular acceptance of the apparatus for the 
PP — ► {pp}s 7T° reaction at 353 MeV and shows a clean n° 
signal with an almost negligible background. Simulations 
indicate that the cm. angular resolution is better than 5°. 

The polarization asymmetry is defined by 



N t /L t - NJL± 



(3) 



v t /x t + Ayv 

where and N± are the numbers of pp — > {pp} s ir° events 
with beam proton spin up and down, corrected for dead 
time, and and L± are the corresponding luminosities. 
The relative luminosity L^/L^ w 0.985 ± 0.015 was esti- 
mated using events at very small polar angles, where the 
polarization asymmetry should be negligible. This proce- 
dure adds about a 3% systematic error to the values of e. 

The analysing power A y is connected to the asymmetry 
through: 



A = 

V P (COS (ppp)' 



(4) 



where P is the transverse polarization of the beam and 
(cos 4>p P ) the average over the diproton azimuthal angular 
distribution. Since the cos 4> pp acceptance is concentrated 
near 1, all the events in the regions analysed contribute 
usefully to the A y measurement. 

The polarization of the proton beam was flipped between 
"spin- up" to "spin-down" (perpendicular to the plane of the 
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Fig. 2. Product of the beam polarisation and analysing power at 
a beam energy of 353 MeV for (a) elastic pp scattering and (b) 
the pp — > dir + reaction. The predictions of the SAID program [21] 
have been scaled to agree with the experimental data and these give 
average COSY proton beam polarisations of (a) P = 0.687 ± 0.008 
and (b) P = 0.668 ± 0.016. In neither case was the uncertainty in 
the SAID prediction included. 



accelerator) every six minutes and no measurements were 
made with an unpolarized beam. The value of P was esti- 
mated from proton-proton elastic scattering and the pp — > 
dit + reaction that were measured in parallel. The analysing 
powers for these reactions were taken from the SAID analy- 
sis program, solutions SP07 forpp —¥ pp and SP96 for pp — > 
dir + [21]. The results of the two methods shown in Fig. 2 
agreed within measurement errors and gave an average po- 
larization of P — 0.68 ± 0.03, where the error includes the 
uncertainties arising from the calibration reactions. 

A simulation was undertaken of the two-dimensional ac- 
ceptance in terms of the pp excitation energy E pp and its 
cm. polar angle 9 pp . This took into account the geome- 
try of the setup and the sensitive areas of the detectors, 
the efficiency of the multiwire proportional chambers and 
the track reconstruction algorithm. In order to avoid po- 
tential problems arising near the limits of the acceptance, 
cuts were made around the edges of the exit window of the 
spectrometer magnet in both the experimental data and 
simulation. This is only a challenge at the larger angles, 
80° < 9^ < 100°, where a compromise had to be made re- 
garding the acceptance ambiguities and this introduces an 
extra 4% systematic uncertainty in this angular region. 

The numbers of detected it events were then corrected 
on an event-by-event basis for acceptance, dead time and 
relative luminosity L^/L^. The latter were important be- 



cause, in the absence of data with an unpolarised beam, an 
average has to be evaluated. 

The luminosity in the experiment was estimated from 
measurements of pp elastic scattering carried out in par- 
allel. The numbers of detected events, corrected for the 
dead time, were compared with a simulation that used a 
generator which included the differential cross section ob- 
tained from the SAID analysis program [21]. Although this 
program does not furnish error bars, experimental data at 
nearby energies suggests that the associated uncertainty is 
about 2%, to which must be added 3% arising from ac- 
ceptance and similar systematic effects. At this level the 
statistical error is negligible and the resulting total lumi- 
nosity was estimated to be 544 ± 22 nb _1 . At this energy 
the pp — > dir + cross section data are less precise than those 
oi pp elastic scattering but, on the basis of the SAID pre- 
dictions, one obtains the completely consistent luminosity 
estimate of 547 nb _1 . 
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Fig. 3. Differential cross section for the pp —y {pp} s vr° reaction at 
353 MeV as a function of the cosine of the pion cm. angle. The solid 
(black) circles represent the ANKE measurements. The errors shown 
here are statistical together with a 4% systematic contribution in the 
80° < 9tt < 100° region coming from the acceptance ambiguity dis- 
cussed in the text. The overall systematic uncertainty is about 4%. 
Open (red) circles are CELSIUS data obtained at 360 MeV [15]. It 
should be noted that the latter data represent averages of measure- 
ments taken in both hemispheres. The curve is a linear fit in cos 2 9^ 
to our data. 

The differential cross section results are presented in 
Fig. 3, where they are compared to those obtained at 
360 MeV at CELSIUS [15]. Within the 10% luminosity 
uncertainty in these data, the overall agreement is very 
good. However, the CELSIUS data at this energy level off 
a little around 90°. This seems to be a feature only of the 
360 MeV results since, at the other energies, linear fits in 
cos 2 On all have good values of x 2 /NDF [15]. 

Fitting our data with a polynomial in cos 2 ^, as in 
Eq. (1), gives parameters 

a = 4.05 ± 0.08 jub/sr, 

a 2 = -2.31 ± 0.14 ,ub/sr, (5) 

Apart from the acceptance uncertainties at the larger 
angles, the error bars quoted here are purely statisti- 
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cal; the ±4% systematic uncertainty from the luminosity 
and acceptance largely cancels in the ratio 02/00. Since 
X 2 /NDF = 23/20, there is clearly no compelling evidence 
for any cos 4 6^ dependence, i.e., a non-zero 04 coefficient, 
and this possibility has been omitted from the curve in 
Fig. 3. 
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Fig. 4. (a) The product of the measured analysing power and differ- 
ential cross section for the pp — ¥ {pp} s tt reaction, with cross section 
errors as in Fig. 3. The curve represents the best fit of Eq. (2), with 
&2 = 1-82 /ib/sr and all higher terms eliminated, (b) Measured values 
of A y for the pp — > {pp} s T reaction. The errors shown are purely 
statistical; the overall systematic uncertainty is about 5%. The line 
represents the quotient of the best fit in pancl-a and the fit to the 
cross section in Fig 3. 



The results for the analysing power of the pp — > {pp} s tt° 
reaction are displayed in Fig. 4, with A y {d<j /dCt) being 
shown in panel a and A y in panel b. These observables must 
be antisymmetric about 90° and the crossing of the data 
through zero around this angle is some confirmation of our 
estimation of L^jL^. These data are subject to the over- 
all uncertainties associated with the luminosity and accep- 
tance evaluation, though these are not relevant for the A y 
in Fig. 4b. There remains, however, the ±3% arising from 
the uncertainty in the value of L^/L^. 

The Ayda/dCl data are consistent with a sin ^ cos #„■ 
behaviour and a fit using the general form of Eq. (2) yields 



b 2 = 1.82 ±0.10 /xb/sr 



(6) 



with x 2 /NDF = 15/21. There is therefore no evidence for 
any sin 6^ cos 3 6^ dependence, i.e., a non-zero 64 coefficient. 
The resulting fits are shown in Fig. 4. 

In order to understand the significance of the results re- 
ported here, we must attempt a partial wave description of 



the data. The most general form of the reaction amplitude 
is 

M = AS ■ p + BS ■ k, (7) 

where S is the polarisation vector of the initial pp spin- 
triplet state, p and k are unit vectors in the cm. frame 
along the directions of the incident proton and final pion, 
respectively. 

The observables studied here are expressed in terms of 
the two scalar amplitudes A and B through [14] 



= ^(\A\ 2 + \B\ 2 + 2Re[AB*}cos8 v ), 



(8) 



The experimental data show no evidence for high par- 
tial waves at 353 MeV and so we model these results with 
only I = and I = 2 contributions. The latter can arise 
from initial L = 1 or L = 3 waves so that, in total, there 
are three possible transitions, 3 P — > 1 Sos, 3 ^2 — > 1 Sod, 
and 3 F 2 — > 1 Sod, see e.g. Ref. [22] for the explicit form of 
the spin-angular structures. We denote the corresponding 
amplitudes by Mf , , and , respectively. 

Expanding the scalar amplitudes in terms of these partial 
waves gives 
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(9) 



Equations (8) and (9) then allow one to relate the mea- 
sured observables of Eqs. (1) and (2) to the partial wave 
amplitudes. For consistency, since we have neglected any 
possible effects arising from s-g interference, we shall also 
drop terms that are bilinear in rf-wave production ampli- 
tudes. In this approximation 



ao = \Mf | 2 - §Re [Mf (Mj 3 + §Mj)*] 
a 2 = 2Re[Mf(M ( f + |M ( f)*], 
fe 2 = 2Im [M?(M?-$M?r], 



(10) 



and so the data only provide three relations between the 
three complex amplitudes. The transverse spin correlation 
parameters contain no extra information since A y>y = 1 
and this is also true for A XtX up to d-d interference terms. If 
the longitudinal-transverse spin correlation parameter A x . z 
were measured, this would provide one further relation but 
this would still not be sufficient for an unambiguous partial 
wave decomposition. For this we need information about 
the phases of the production amplitudes. 

The 3 P partial wave is uncoupled and, at the energy 
where the experiment was performed, its inelasticity is very 
small. Under these conditions the Watson theorem, which 
fixes the phase induced by the initial state interaction to 
that of the clastic proton-proton scattering, applies [23]. 
Thus we take M s p = |Mf \e i53p o , with 6s Po = -14.8° [21]. 
Note that we do not include any phase associated with the 



A 



1 Sq final pp state because it is common for all partial waves 
and therefore does not affect the observables. 

For coupled channels, such as 3 P 2 — 3 F 2 , the strict condi- 
tions of the Watson theorem do not apply. However, at our 
energy the mixing parameter, as well as the inelasticities, 
are still negligibly small and thus to a good approximation 
we may also use the Watson theorem here. Further evidence 
for the smallness of the channel coupling is to be found in 
two potential models [24,25]. In both models the T-matrix 
for the transition from the 3 F 2 to the 3 P 2 wave is almost 
real; the phase of M/j* is driven by 5-ip 2 — 17.9°, whereas 
the phase of Ml can be neglected. The quality of this ap- 
proximation was also checked by explicit calculations of the 
<i-wave production amplitudes within chiral effective field 
theory. These showed that up to order m^/inN the above 
phase assumptions should be valid to within ±2°. 

Using the phase information in this way, we find that 

Mf= (55.3 ± 0.4) - (14.7 ± 0.1)i ^nb/sr , 

Mf = -(26.6 ±1.1) - (8.6±0.4)i ^nb/sr, 

Ml = (5.3 ±2.3) Vnb/sr. (11) 

The values quoted here were obtained by considering also 
our np — > {pp} s n~ data though the numbers would change 
but marginally if one included only the pp — > {pp} s n° re- 
sults in the fit. The error bars quoted here are statistical and 
do not include the overall systematic uncertainties. How- 
ever, changing the normalisations of the differential cross 
section and analysing powers by 3% and 4%, respectively, 
leads to changes that are comparable to the quoted errors. 
On the other hand, we could not investigate the less tan- 
gible ones associated with the neglect of the channel cou- 
pling and the truncation in the partial wave expansion. The 
weakness of pion production from the initial 3 F 2 waves at 
353 MeV, in addition to being in agreement with theoret- 
ical prejudices, is also consistent with the low inelasticity 
found for this wave [21]. 

In summary, we have measured the differential cross sec- 
tion and analysing power of the pp — > {pp} s ir° reaction at 
353 MeV. The angular distributions of A y and da/dSl are 
both well represented by retaining only pion s and d waves 
in a phenomenological description. The values of da/dQ 
agree well with the results obtained at CELSIUS [15] over 
most of the angular range. However, at this energy these 
data flatten off around the middle of the angular distribu- 
tion and, if this effect were correct, it would signal a large 
contribution from 3 F 2 or even higher partial waves. 

By making plausible assumptions on the coupling be- 
tween the nucleon-nucleon channels and invoking the Wat- 
son theorem it was possible to estimate the partial wave 
amplitudes with their phases. These could be checked 
through a future measurement of the spin-correlation pa- 
rameter A XfZ , though this would require the installation of 
a Siberian snake in COSY. 

In an associated letter [26] , the isospin-0 amplitudes are 
investigated through the measurements of the differential 
cross section and analysing power of the quasi-free pn — > 



{pp} s ir~ reaction in this energy domain. The extraction of 
the p-wave amplitudes from the data of Ref. [26] required 
a knowledge of the s- and d-w&ve amplitudes of the type 
provided here. In addition, data have already been taken 
on the transverse spin correlation parameter A XiX for this 
reaction [27]. The full collection of these results will lead 
to very useful constraints on the parameters of the chiral 
effective field theory that link pion production to the three- 
nucleon interaction [1,2]. 
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